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Recently, a phase transition to synchronized congested traffic has been observed in empirical highway 
data [B. S. Kerner and H. Rehborn, Phys. Rev. Lett. 79, 4030 (1997)]. This hysteretic transition 
has been described by a non-local, gas-kinetic-based traffic model [D. Helbing and M. Treiber, Phys. 
Rev. Lett. 81, 3042 (1998)] that, however, did not display the wide scattering of synchronized 
states. Here, it is shown that the latter can be reproduced by a mixture of different vehicle types 
like cars and trucks. The simulation results are in good agreement with Dutch highway data. 
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Recent publications stressed the fact that, whereas in 
free traffic the observed flow-density diagram is well de- 
scribed by a unique one-dimensional flow-density rela- 
tion, in congested traffic the empirical data points are 
rather distributed over a two-dimensional region QQ] . In 
order to account for this fact, Kraufi || has recently pro- 
posed that the driver behavior is changed in congested 
traffic compared to free traffic. Instead of this, Lenz et 
al. J|| and Schreckenberg (5) have suggested that the wide 
scattering of data is caused by an anticipation effect of 
drivers who not only react to the respective vehicle in 
front but also to the traffic dynamics further ahead. 

In contrast to these "microscopic" approaches which 
simulate the interactions of individual vehicles, macro- 
scopic models describe the evolution of the macroscopic 
velocity V[x, t) = (v a ) and the vehicle density p{x,t) = 
(1/sa), which are local averages of the "microscopic" ve- 
locities v a of the vehicles a and their center-to-center 
distances s a . All fluctuating quantities like individual ve- 
locity variations or distance distributions are eliminated. 
This means that, in deterministic macro-simulations, all 
self-organized structures (like stop-and-go waves or con- 
gested traffic) are smooth. 

Therefore, some researchers believe that, while the 
wide scattering of the congested flow-density data may 
be reproduced by microsopic traffic models, macroscopic 
ones will fail for principal reasons. However, motivated 
by the circumstance that the scattering has been ob- 
served in aggregated rather than single-vehicle data, we 
are confident that a macroscopic simulation of this effect 
should be possible. 

In the following, we will show that the scattering can be 
explained by the fluctuations caused by a heterogeneous 
traffic population, which enter the macroscopic simula- 
tions via the boundary conditions. We will distinguish 
cars and trucks characterized by different sets of param- 
eter values. These define two equilibrium flow-density 
relations of pure car traffic and pure truck traffic, re- 
spectively, which are close to each other at small vehicle 
densities, but considerably different in the congested den- 
sity regime. For mixed traffic, we interpolate between 
both parameter sets and, hence, between both equilib- 
rium relations, using a weighted average. The weights 



are extracted from real traffic data by determining the 
proportion of long vehicles ("trucks"). This method al- 
lows to simulate a uni-directional multi-lane freeway by 
an effective one-lane model for one car species with av- 
erage, but varying parameter values. Our simulations 
are carried out with the empirically measured boundary 
conditions, and the results are quite realistic. 

For the simulations, we use the macroscopic, gas- 
kinetic-based traffic model (GKT model) [glJj, which 
shows a realistic instability diagram and the characteris- 
tic properties of traffic flows demanded by Kerner and 
Rehborn |IJ] . More importantly, this model is able to de- 
scribe the hysteretic phase transition to congested states 
with high traffic flow || (called "synchronized traffic" 
Q) which typically occurs behind on-ramps, gradients, 
or other bottlenecks of busy freeways j(| . 

According to the gas-kinetic-based model, the evolu- 
tion of the vehicle density p(x, t) in dependence of the 
time t and the position x along the freeway is given by 
the continuity equation 
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Here, V(x, t) denotes the average velocity of the vehicles. 
At on-ramps (or off-ramps), the source term Q r mp/{nL) 
is given by the actually observed inflow Q rm p > from 
(or outflow Q rmp < to) the ramp, divided by the merg- 
ing length L and by the number n of lanes. Otherwise 
it is zero, reflecting the conservation of the number of 
vehicles. The average velocity obeys the equation 



dV 



v 



dV 

dx 



1 d( P 9) 
p dx 



1 



TransportTerm p rC ssurcTerr 



+ -{U-V) 



RelaxationTcrm 



(2) 



According to this, the temporal change of the average 
velocity is given by a transport term (caused by a propa- 
gation of the velocity profile with V), a so-called pressure 
term (that reflects dispersion effects due to the finite ve- 
locity variance 9 of the vehicles) , and a relaxation term 
(describing the adaptation to a dynamic equilibrium ve- 
locity U with a certain relaxation time r). In our gas- 
kinetic-based model, the analytically derived formula for 
the dynamical equilibrium velocity is 
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where Vq is the desired (maximum) velocity, T the aver- 
age time headway at large densities, and p max the max- 
imum vehicle density. A prime indicates that the corre- 
sponding variable is taken at the advanced "interaction 
point" x' = x+j(l/ p max +TV) rather than at the actual 
position x. This accounts for the anisotropic anticipa- 
tion behavior of drivers. The monotonically increasing 
"Boltzmann factor" 
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can be derived from gas-kinetic formulas and describes 
the dependence of the braking interaction on the dimen- 
sionless velocity difference 5y = (V — V')/y/9 + 9' . Fi- 
nally, the dynamics of the variance can be approximated 
by the constitutive relation 
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where the coefficients A = 0.008, AA = 0.015, p c = 
0.28p max , and Ap = 0.1p max have been obtained from 
single- vehicle data @ . 

The velocity-density relation resulting for this model in 
spatially homogeneous and stationary equilibrium reads 



V 2 

™ = 2V Q 



with 
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This also determines the equilibrium traffic flow by 
Qe(p) = pV e (p) , 



(7) 
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which, for a given parameter set, is a one-dimensional 
curve. However, as will be shown in the following, the 
empirically observed two-dimensional region of " synchro- 
nized" congested states can be reproduced by simulating 
a mixture of different vehicle types. Although it has 
not been stressed clear enough, it is known from mi- 
crosimulations that heterogeneous traffic produces con- 
siderable fluctuations of the aggregate quantities like the 
vehicle density and the average velocity |]-[ll|]. Never- 
theless, we do not need to carry out microsimulations to 
account for the two-dimensional scattering of synchro- 
nized traffic states. It is sufficient to simulate traffic in 
a macroscopic way with empirically obtained boundary 
conditions, including the varying proportion of long ve- 
hicles ("trucks"). Thus, we do not need to assume other 



sources of fluctuations than observed ones. A reasonable 
agreement with empirical data can already be reached 
by distinguishing two vehicles types only, short vehicles 
(" cars" ) and long ones (" trucks" , with a length of at least 
7 m). Each type is characterized by its own paramter set. 
For the cars we assume a desired velocity Vb = 112km/h, 
an average time headway T = 1.0 s at large densities, and 
a maximum density p max = 110 vehicles/km. Trucks are 
described by the parameters Vq = 90km/h, T — 5 s, and 
Pmax = 100 vehicles/km. The remaining model parame- 
ters are the same for both types: r = 25 s and 7 = 1.6. 
The parameters in the constitutive relation (g|) for the 
variance have also been chosen identical. 

According to the philosophy of macroscopic models, we 
now define time-dependent " effective" model parameters 
X{t) as weighted averages of the respective car and truck 
parameters A car and X tluc k' 



X (t) = Ptruck(t)X tTuck + [1 - Ptruck(t)]X c 



(9) 



Here, ptruck(i) is the proportion of trucks averaged over 
a time interval At around t [Figure |(a)] . Although the 
approximations behind the resulting "effective" macro- 
scopic simulation model are rather crude, it yields a sur- 
prisingly good agreement with empirical data. Even bet- 
ter results are expected for macroscopic models which 
explicitly take into account different vehicle types and 
lane-changing interactions among the freeway lanes flT2| . 

We simulated traffic flow on a section of the Dutch two- 
lane motorway A9 from Haarlem to Amsterdam (Fig- 
ure |lj) from the detector cross-section Dl (0 km) to D6 
(5.7km). For this purpose, the measured single-vehicle 
data were aggregated to 1-minute averages of the veloc- 
ity, traffic flow, and truck proportion (Figure ^) . Between 
7:30 am (450 min) and 9:30 am (570 min) in the morn- 
ing of November 2, 1994, we find transitions from a low- 
density regime to a high-density regime corresponding to 
transitions between free and congested traffic. Figure || 
illustrates that the congested state at D2 is connected 
with a considerable velocity drop, while the flow is de- 
creased only by about 10%, both in the empirical data 
and in the simulation. In addition, the congested traffic 
state relaxes to free traffic downstream of the on-ramps 
[Figures f|(c) and (d)]. A comparison with Figures 1(b) 
and 3(c) of Ref. M suggests that the congestion in the 
investigated data corresponds to synchronized traffic. 

As inflow and outflow boundary conditions, we used 
the data of the cross-sections Dl and D6, respectively, as 
shown in Figure |^(b). There are two on-ramps and one 
off-ramp in the considered section. For all ramps, we use 
the empirical data of the traffic flow Q rmp , divided by 
the number n — 2 of lanes |(| [Figure ||(c)] , and assume 
a merging length L = 200 m. 

In the simulation, congested traffic first sets in at 
t w 450 min near the on-ramp at D3a, which agrees well 
with the empirical findings. We started the simulation 
50 min earlier to eliminate any effects of initial condi- 
tions, and to show the spontaneous nature of the tran- 
sition [Figure ||(a)] . In Figure ^(b), the free traffic flow 
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before the breakdown (t < 450 min) and after the re- 
covery (t > 570 min) is delineated by the points at the 
low-density (left) side of the diagram, which more or less 
define a one-dimensional curve. In contrast, the con- 
gested traffic state is represented by the points at the 
high-density (right) side, which are distributed over a 
two-dimensional region. Some minutes later, the front 
of the congested state crosses the on-ramp at D2, which 
causes congested traffic upstream of it. In accordance 
with the mechanism of the formation of synchronized 
traffic proposed in the congested state upstream of 
D2 has a lower flow and a higher density [Figure |](a)] 
than that between D2 and D3a. The congested states 
are sustained for nearly two hours, until the inflows from 
both the main road and the on-ramps are considerably 
decreased, which shows the hysteretic nature of the tran- 
sition. 

Summarizing our results, one can say that the macro- 
scopic, gas-kinetic-based traffic model allows to simulate 
synchronized traffic, including the associated scattering 
of the flow-density data in the congested regime. Sim- 
ulations of this model with only one vehicle type PJTsJ 
suggested that the phenomenon of synchronized traffic 
as such (i.e., high traffic flows at low velocities) does 
not depend on the existence of different types of vehi- 
cles. However, as is often the case for self-organized non- 
chaotic patterns resulting from deterministic dynamics, 
the flow-density diagram is essentially one-dimensional. 

In this Letter, we showed that a realistic scattering in 
the flow-density plane can be simulated by distinguishing 
several vehicle types with different parameter sets, the 
measured proportions of which are the weights for deter- 
mining the time-dependent "effective" parameter set. A 
reasonable agreement with empirical data from Dutch 
highways is already obtained for two different vehicle 
types, cars and trucks. Our results also indicate that, 
when studying dynamical phenomena in empirical traf- 
fic data, it is highly recommended to thoroughly analyze 
the proportion of trucks, which shows surprisingly large 
variations [Figure §(a)]. 

Notice that the assumed parameter variations due to a 
changing truck fraction can explain both the relatively 
low scattering of flow-density data in the low-density 
regime and the wide scattering in the regime of congested 
traffic. While the small amount of scattering for free traf- 
fic at low densities is caused mainly by variations of the 
individual desired velocity (with a standard deviation of 
about 10% of the mean value), the main reason for the 
considerable scattering for congested traffic at densities 
above 30 vehicles/km are the variations of the time head- 
way (which are of the order of 100%). There are other 
effects that influence scattering of congested traffic (e.g., 
lane changes), but they will only increase the scattering 
of the flow-density data. 

Finally, we mention that the equations of the gas- 
kinetic-based traffic model together with relation (^) 
for the stochastic quantities V c , T, and p max , represent 
stochastic partial equations with multiplicative noise. 



They may serve as prototype for introducing stochastic- 
ity into macroscopic equations in a controlled and empir- 
ically justified manner. 
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FIG. 1. Overview of the evaluated stretch of the Dutch 
Highway A9 from Haarlem to Amsterdam. 
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FIG. 2. (a) Proportion of trucks, from one-minute av- 
erages, (b) Upstream and downstream boundary conditions 
for the flow, taken from measured one-minute data at the 
cross-sections Dl and D6, respectively, (c) Flows of the three 
ramps in the considered section. The off-ramp at detector 
Dl was left out. It leads only to changes of the traffic situa- 
tion upstream of the cross-section Dl for which no data were 
available. 
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FIG. 3. (a) Velocity, and (b) traffic flow at D2 according to 
the model, in comparison with the empirical one-minute data. 
The breakdown of velocity is a result of a dynamical transi- 
tion, since neither the initial conditions, nor the boundary 
conditions, or the ramp flows used in the simulations contain 
any significant peaks. 
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FIG. 4. The displayed points in density-flow space corre- 
spond to empirical 1-minute data (dark crosses) and related 
simulation results (grey boxes), separately for the cross sec- 
tions D2, D3, D4, and D5. The simulations manage to re- 
produce both the quasi-linear flow-density relation at small 
densities and the scattering over a two-dimensional region at 
high densities. For comparison, we have displayed the equilib- 
rium flow-density relations for traffic consisting of 100% cars 
(— ), and 100% trucks ( ). 
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